Optical absorption spectra in SrCi^Os two-leg spin ladder 
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We calculate the phonon-assisted optical-absorption spectra in SrCu2C>3 two-leg spin- ladder sys- 
tems. The results for two models proposed for SrCu2C>3 are compared. In the model including the 
effects of a cyclic four-spin interaction, the shoulder structure appears at ~ 978cm -1 and the peak 
appears at ~ 1975cm -1 in the spectrum for polarization of the electric field parallel to the legs. 
In the other model which describes a pure two-leg ladder, the peak appears around the lower edge 



of the spectrum at 
SrCu 2 3 . 



1344cm" 1 . The feature can be effective in determining the proper model for 



PACS numbers: 75.40.Gb, 78.30.-j, 74.72. Jt, 75.40.Mg 

Intensive studies for 5 = 1/2 two-leg spin-ladder sys- 
tems have revealed fascinating aspects of elementary ex- 
citations as well as thermodynamic properties of the sys- 
tem p|. Using several theoretical methods, it was shown 
that 5 = and 5 = 1 two-triplet bound states lie below 
the two-triplet continuum in addition to the one-triplet 

mode HHSEH0EI!,[l3. Furthermore : the effects 
of a cyclic four-spin interaction on the 5=1/2 two- leg 
spin-ladder systems have attracted a great amount of at- 
tention recently [HI E3, Q [IE llE E3, EE E3, 
l2lj . From the analysis of the experimental results for 
La6CasCu2404i observed by inelastic neutron-scattering 
experiments, it was corroborated that a cyclic four-spin 
interaction is necessary to explain the observed disper- 
sion relation of the one-triplet mode . 

For SrCu203, which is considered to be a typical 
5 = 1/2 two- leg spin-ladder material, two models were 
proposed to reproduce the temperature dependence of 
the susceptibility. One model includes the effects of a 
cyclic four-spin interaction , while the other model de- 
scribes a pure two-leg ladder [22J. It seems, thus, difficult 
to decide the proper minimal model only from thermo- 
dynamic quantities such as the susceptibility. On the ba- 
sis of the observations, dynamical structure factors have 
been investigated for the two models 0, . Although 
the dispersion relation for the model including the ef- 
fects of a cyclic four-spin interaction becomes flatter, the 
difference in the distribution of the weights between the 
two models may be difficult to be observed by inelas- 
tic neutron-scattering experiments [l^. Other methods 
to detect characteristics of dynamical properties are re- 
quired to decide the proper model for SrCu203. 

The midinfrared optical-absorption spectra measured, 
e.g., by optical conductivity in low-dimensional cuprates 
can be successfully analyzed in terms of the phonon- 
assisted absorption mechanism [23-31, 9, 10]. In this 
mechanism, an absorbed photon simultaneously creates 
a Cu — O bond-stretching phonon and a multitriplet ex- 
citation with a singlet coupling. This simultaneous ex- 
citation enables us to probe a phonon-weighted sum of 
the 5 = multitriplet excitation from the entire Bril- 
louin zone. In fact, it was confirmed from the optical- 
conductivity measurements that the 5 = two-triplet 



bound state exists below the two-triplet continuum in a 
5 = 1/2 two-leg spin-ladder material (Ca, La)i4Cu2404i 
[31, 9, 10]. 

In this paper, we investigate the optical-absorption 
spectra caused by the phonon-assisted two-triplet exci- 
tation in SrCu203 two-leg spin-ladder systems on the 
basis of the proposed two minimal models. Using a con- 
tinued fraction method based on the Lanczos algorithm 
[23, we first calculate the wave- number-resolved 5 = 
two-triplet excitation spectra . We then calculate the 
optical-absorption spectra and discuss characteristics of 
the results for the two models of SrCu203. 

We start our discussion with the following Hamilto- 
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N/2 

+ Jcyc y"^{4[QSl,i ' 5i i i+i)(52,i+i ' ^2, i) 
t=l 

+ {Sl,i • S2,i)(,Si l i+i ■ S>2,i+l) 
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N/2 

+ Jdiag ' ^2,t+l) + (5l,i+l ' ^2,-t)]; (1) 

i=l 

where 5^, denotes the 5 = 1/2 spin operator in the ith 
rung of the I = 1,2 chain, N is the total number of 
spins, and J|| and Jj_ are the coupling constants along the 
leg and rung, respectively. The coupling constants for a 
cyclic four-spin interaction and a diagonal interaction are 
denoted as J cyc and Jdiag) respectively. The coupling con- 
stants used in the two models for SrCu203 are (i) J± — 



150meV, Ji, = 195meV, J, 



eye 



18meV, Jdia g = 3meV 



(Ref. 11) and (ii) J ± = 86mcV, J\\ = 172meV, J cyc = 
^diag = (Ref. 22). The periodic boundary condition is 
applied along the chain. 
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The wave-number-resolved 5 = two-triplet excita- 
tion spectra for polarization of the electric field parallel 
to the legs (E\\ legs) can be expressed as 



R l ^{q x ,q y ;uj) 



1 -d'eg 

z-n 



(2) 



where (fc and %(= 0, 7r) are the wave numbers along the 
leg and rung, respectively, l^o) is the eigenfunction of 
the lowest eigenvalue Eq, and z = u> + is + Eq. The 
operator A 1 ^ = {l/^^ .e^+^S^ ■ S l<j+1 is 
the Fourier transform of the operator for the locally ex- 
citing 5 = two-triplet state on the neighboring sites 
in a chain. We set fi = 1 and e — 3.0 x 10~ 2 . The 
energy is measured in units of Ju . In the expression 
(2), C le6 (q x ,q y ;uj) can be represented in the form of 
the continued fraction, which can be calculated numer- 
ically by a Lanczos algorithm [33|. The total contribu- 
tion of C leg (q x , q y ; uj) for fixed q x and q y is normalized to 
unity, because the following sum rule has to be satisfied: 
F lee (qx,q y ) = Jq 00 dw-R les ((fc, q y ;uj). We use the two sets 
of coupling constants corresponding to the two models 
for SrCu 2 3 : (i) JJJ\\ = 0.769, J cyc / J\\ = 0.0923, and 
J diag /J|l = 0.0153 (Ref. 11) and (ii) J±/J\\ = 0.5 and 

Jcyc = Jdi&g = (Ref. 22). 

In Fig. 1, R les (q x , 0;u) and R lcg (q x , tt; uj) in the two 
models are shown for N — 28. The weight is proportional 
to the area of the full circle. In model (i) the convergences 
have the relative errors of about O(10~ 10 ) for uj < 3 and 
O(10~ 2 ) - 0(1O~ 5 ) for ui > 4, while in model (ii) they 
have the relative errors of about O(10~ 12 ) for uj < 3 and 
O(10~ 3 ) - O(10~ 5 ) for oj > 4. In both models, the large 
weights for given q x lie in the lowest excited states in 
q x > 0.4tt. 
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FIG. 1: (a) R lcg (q x ,0;uj) and (b) R leg -(q x , tt; uj) in model (i) 
and model (ii) for N = 28. The weight is proportional to the 
area of the full circle. 



To discuss whether such lowest excited states form a 
lower edge of the excitation continuum or an isolated 



mode of the bound state, we next investigate the finite- 
size effects of the poles and their residues (33 of the con- 
tinued fraction C lcg (q x , q y ;uj) 0, As reported in 
Refs. 34-36 and 18, a pole which belongs to an excitation 
continuum tends to have an appreciable size dependence 
of at least either its position or its residue. Since the size 
dependence of the position of the pole is hardly seen, we 
only show the size dependence of the residue in the fol- 
lowing. In Fig. 2, we first show R lcs (q x , 0; uj) and the size 
dependence of the lowest excited states of C les (q x ,0;uj) 



for J±/J\\ = 1 and J c 



^diag = 0. In q x > 0.57T, the 
size dependence is scarcely seen. Thus, the lowest ex- 
cited states in q x > 0.57T form an isolated mode of the 
5 = two-triplet bound state, while those in q x < 0.5n 
form the lower edge of the excitation continuum. The 
results are consistent with those obtained for the other 
two methods H, . 
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FIG. 2: R lcs (q x , 0; u>) for TV = 28 and the size dependence of 
the residues of the lowest excited states for C le£ (q x , 0; u>) in 



J±/J\\ = 1 and J c - 



Jdia 



: 0. 



In Fig. 3, we show the size dependence of the residue 
for the lowest excited states. In C lcg (q x , tt; uj), the 
residues for both models have a noticeable size depen- 
dence in < q x < tt, indicating that the lowest excited 
states of the q y = it mode become the lower edges of the 
excitation continuum in the thermodynamic limit. In 
C lcs (q x , 0; uj) for model (ii), on the other hand, the size 
dependence in q x > 0.57T is very weak, while the residue 
in q x < 0.57T decreases with increasing N. Therefore, 
the lowest excited states of the q y = mode in model 
(ii) form an isolated mode in q x > 0.5n, and those in 
q x < 0.57T become the lower edge of the continuum. In 
C les (q x , 0; uj) for model (i), the residues around q x ~ 0.67T 
and 0.87T decrease slightly with increasing N. Although it 
is difficult to draw a definite conclusion, the observations 
may indicate that the lowest excited states of the q y = 
mode in model (i) become the lower edge of the excita- 
tion continuum in q x > 0.67T. In fact, the size dependence 
of the residue in q x > 0.6n is quite similar to that for the 
5 = 1/2 isotropic Heisenberg chain in q > 0.67T [35| . 
where the lowest excited states form the lower edge of 
the two-spinon continuum. Note that in Ref. 10, it is 
argued that the 5 = two-triplet bound state exists in a 
larger J cyc . To develop more definite finite-size analysis, 
a larger system has to be calculated. 

The optical-absorption coefficient for E || legs is given 
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FIG. 3: The finite-size effects of the residues of the lowest 
excited states for (a) C lcg (q x , 0; lo) and (b) C les {q x ,n;ui) in 
model (i) and model (ii). 
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by HElllllli 

a les (u>) = a o L0l les (u> - w ), (3) 

where a is a constant depending on the material, lo q is 
the frequency of the Cu — O bond-stretching phonon, and 
^ log M = E g ^4^ los (fe,9<,;w) with the phonon 
form factors = sin 4 (^2) and = sin 2 (q x /2) + 
sin 4 (g x /2). To calculate P CK (u — u ), we apply the spline 
interpolation to the lowest excited states of R l ° s (q x , q y \Lo) 
in the q x direction. 

In Fig. 4, we show I 1 ° s (lo — u> ) thus obtained for the 
two models in N — 28. The dashed line and dotted 
line represent the spectra from the q y — mode and 
q y = 7r mode, respectively, and the solid line indicates 
their sum. A noticeable difference between the results 
for the two models appears. In I lcg (ui — uj ) of model 
(i), the shoulder appears at u> — lo q ~ 0.7 and the peak 
appears at u - w ~ 1.2. In I lcg (uj — uj ) of model (ii), 
on the other hand, the peak appears at the lower edge of 
the spectrum: u> — u> ~ 1.0. In model (i), the dispersion 
relation for the lowest excited states of the q y = mode 
becomes flatter, which yields the peak at w - w ~ 1.2. 
The shoulder structure is caused by the lowest excited 
states of the q y — ir mode, which lies below those of 
the q y = mode in q x > 0.87T. Note that when we set 
Jcyc — in model (i), the lowest excited states of the 
q y = 7r mode lie above those of the q y = mode and no 
shoulder structure emerges. Therefore, such shoulder and 
peak structures in model (i) are caused by the cyclic four- 
spin interaction. In model (ii), the dispersion relation of 
the lowest excited states in the q y = mode lies below 
that of the q y = ir mode. Thus, in model (ii) no shoulder- 
appears and the peak structure is seen at the lower edge 
of the spectrum. 

In Fig. 5(a) we show the size dependence of the ex- 
citation energies for the shoulder and peak, and in Fig. 
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FIG. 4: The optical-absorption spectra for E \\ legs in model 
(i) and model (ii) for TV = 28. The dashed lines and dotted 
lines represent the spectra from the q y = mode and q y = n 
mode, respectively, and the solid lines indicate their sum. 



5(b) we show the size dependence of their weights. As 
mentioned above, the shoulder is caused by the q y = ir 
mode (see the dotted lines in Fig. 4) and the peak is 
caused by the q y = mode (see the dashed lines in Fig. 
4). In model (i), the extrapolated excitation energies for 
the shoulder and peak at N — > oo are lo — lo ~ 0.240 and 
lo — lu ~ 0.874, respectively. The extrapolated weights 
of the shoulder and peak at N — > oo are ~ 0.0996 and 
~ 0.474, respectively. Therefore, in model (i) the shoul- 
der and peak structures in the optical-absorption spec- 
trum are intrinsic. In model (ii), the extrapolated ex- 
citation energies for the shoulder and peak at N — > oo 
almost coincide and take the values lo — lo ~ 0.540 and 
0.536 with corresponding weights of ~ 0.0553 and 0.876. 
Thus, in model (ii) only one peak appears around the 
lower edge in the optical-absorption spectrum. 

Using J|| = 195meV in model (i) HJ, J|| = 172meV in 
model (ii) l22l . and lo = 600cm" 1 commonly evaluated 
in cuprates j3J, the excitation energies for the shoulder 
and peak in the model (i) can be evaluated as 978cm _1 
and 1975cm , respectively, while the peak in the model 
(ii) can be evaluated as 1344cm -1 . Since a leg (o;) can be 
observed experimentally by the optical conductivity, such 
different characteristics between the optical-absorption 
spectra in the two models are probably detectable. The 
observed findings can be effective in determining the ap- 
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FIG. 5: The size dependence of (a) the excitation energies 
for the peaks (see the dashed lines in Fig. 4) and shoulders 
(see the dotted lines in Fig. 4), and of (b) their weights. 



propriate minimal model for SrCu203. 

In summary, we have investigated the phonon-assisted 
optical-absorption spectra on the basis of the two models 
proposed for SrCu2 03. Characteristics of the lowest ex- 
cited states have been discussed in connection with the 
shoulder and peak structures in the optical-absorption 
spectra. In the model including the effects of a cyclic 
four-spin interaction, the shoulder and peak structures 
appear at ~ 978cm -1 and ~ 1975cm -1 , respectively, 
in the optical- absorption spectrum for E || legs. In the 
other model which describes a pure two-leg ladder, no 
shoulder appears and the peak appears around the lower 
edge of the optical-absorption spectrum for E\\ legs at 
~ 1344cm -1 . 

Our computational programs are based on TITPACK 
version 2 by H. Nishimori. Numerical computations were 
carried out at the Yukawa Institute Computer Facility, 
Kyoto University, and the Supercomputer Center, the In- 
stitute for Solid State Physics, University of Tokyo. This 
work was supported by a Grant-in- Aid for Scientific Re- 
search from the Ministry of Education, Culture, Sports, 
Science, and Technology, Japan. 
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